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We study possible effects of CP violation in the Higgs sector on tt production at a 77- 
collider. These studies are performed in a model-independent way in terms of six form-factors 
{5R(<S' t ), K(P 7 ), 7(P 7 ), St, Pt} which parametrize the CP mixing in the Higgs sector, and a 

strategy for their determination is developed. We observe that the angular distribution of the decay 
lepton from lft produced in this process is independent of any CP violation in the tbW vertex and 
hence best suited for studying CP mixing in the Higgs sector. Analytical expressions are obtained 
for the angular distribution of leptons in the c.m. frame of the two colliding photons for a general 
polarization state of the incoming photons. We construct combined asymmetries in the initial state 
lepton (photon) polarization and the final state lepton charge. They involve CP even (*’s) and odd 
( y’s ) combinations of the mixing parameters. We study limits up to which the values of x and y, 
with only two of them allowed to vary at a time, can be probed by measurements of these asymme¬ 
tries, using circularly polarized photons. We use the numerical values of the asymmetries predicted 
by various models to discriminate among them. We show that this method can be sensitive to the 
loop-induced CP violation in the Higgs sector in the MSSM. 

PACS numbers: 12.60Fr, 14.65Ha, 14.80 Cp, 11.30 Pb. 


I. INTRODUCTION 


The Standard Model (SM) has been tested to an 
extremely high degree of accuracy, reaching its high 
point in the precision measurements at LEP. However, 
the bosonic sector of the SM in not yet complete, the 
Higgs boson is yet to be found. A direct experimental 
demonstration of the Higgs mechanism of the fermion 
mass generation still does not exist. Also lacking is a first 
principle understanding of CP violation in the SM. In 
this note we look at the possibilities of probing potential 
CP violation in the Higgs sector at the proposed 77 
colliders [1|. 

Such a study necessarily means that we are look¬ 
ing at models with an extended Higgs sector. CP 
violation in the Higgs sector can be either explicit, one 
of the first formulations of such a CP violation being the 
Weinberg Model or can be spontaneous ||, where 
the vacuum becomes CP non-invariant. The mechanism 
for creating CP violation in the Higgs sector could be 
different in different models but all such mechanisms will 
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result in CP mixing and then the mass eigenstate scalar 
will have no definite CP transformation property. In 
specific models with an extended Higgs sector, such as 
the Minimal Supersymmetric Standard Model (MSSM), 
for example, the lightest Higgs remains more or less a CP 
eigenstate and the two heavier states H, A, which would 
be CP-even and CP-odd respectively in the absence of 
CP violation and are close in mass to each other, mix. 
The expected mixing can actually be calculated as a 
function of the various parameters of the model §§ El- 
In our study, however, we do not stick to a particular 
model of CP violation and adopt a model-independent 
approach to study the effects of this CP violation on tt 
production in 77 collisions. Such an approach has been 
adopted in earlier studies 0. 

We study 77 —> tt through the diagrams shown in 
Fig. 0 - It has been observed earlier [flj that there exists 
a polarization asymmetry of the tt produced in the final 
state if the scalar (j> exchanged in the s-channel is not a 
CP eigenstate. We parametrize the V 77? j, V t i^> vertices in 
a model-independent way in terms of six form factors to 
include the CP mixing, following Ref. [ij. We investigate 
in this study the effect of such a CP mixing on the 
angular and charge asymmetries for the decay leptons 
coming from the t/t which reflect the top polarization 
asymmetries. It is known Q that 77 colliders can 
provide crucial information on the CP property of 
the scalar produced in the s-channel, due to the very 
striking dependence of the process on the polarization 
of the 7’s. These colliders will also offer the possibility 













2 



to assess the sensitivity of this process to the size of 
various form factors involved in the parametrization of 
V 77 ? j,V(t 0 . At times we have used specific predictions 
for the form factors in the MSSM |j as a guide and for 
purposes of illustration, in our analysis. We show that 
this process is capable of probing the MSSM loop effects 
using these asymmetries. 


FIG. 1: Feynman diagrams contributing to 77 —> tt produc¬ 
tion. 


of measuring the two-photon width of the SM Higgs 
very accurately § The 77 production of a scalar 
followed by its decay into a Z pair is shown to provide 
crucial information required for a model independent 
confirmation of its spin and parity |ll|. Possibilities of 
studying the MSSM Higgs bosons in 77 collisions in the 
bb and neutralino-pair final states, are shown 0 to 
give access to regions of the supersymmetric parameter 
space not accessible at other colliders. Thus in general 
the 77 colliders will provide a very good laboratory 
for studying the scalar sector. Here we concentrate on 
the polarization asymmetries of the final state t and t 
caused by such a CP violation [|7|j. The large mass of t 
implies that it decays before hadronization. As a result 
it acts as a heavy lepton and the spin information gets 
translated to distribution of the decay leptons. Thus 
we can use these angular distributions as a probe for 
possible CP violation. We use only the decay lepton 
angular distributions and construct asymmetries that 
reflect the t/t polarization asymmetries caused by the 
CP violation in the Higgs couplings. We observe that 
these are independent of any CP-violating contribution 
in the tbW vertex. The same is not true of the energy 
distribution of the decay lepton. Hence we restrict our 
analysis to the angular distributions and keep the tbW 
vertex completely general, choosing the tt'j vertex to be 
standard. The latter of course is relevant only for the 
continuum background. 


The plan of the paper is as follows. In section || 
we give the general form for CP-violating vertices 
involving the Higgs, the t quark and the photons as 
well as the tbW decay vertex for the t quark. The tt 
production and f-decay helicity amplitudes obtained 
with these vertices are then presented. In section 

III we obtain an analytic expression for the angular 
distribution of the decay lepton in t/t decay. We discuss 
the insensitivity of the decay-lepton angular distribution 
to the anomalous coupling in the tbW vertex. Section 

IV deals with the ideal photon collider [0. Numerical 
results are presented in sectio n M , discussed in section 
VI and we conclude in section VII. 


II. INTERACTION VERTICES AND 
HELICITY AMPLITUDES 


The interaction vertex of t with a scalar cf >, which may 
or may not be a CP eigenstate, may be written in a 
model-independent way as 

Vut = ~ie^{S t +n 5 P t ). ( 1 ) 


The general expression for the loop-induced 77 <fi vertex 
can be parametrized as 


V 77 <^ — 


—i^fsa 

47T 


S-y(s) ( ei.e 2 -(ei.fc 2 )(e 2 A’i) 


We develop a strategy to study the CP property 
of the Higgs by looking at angular distributions of 
leptons and antileptons, for different polarizations of the 
colliding photons. Towards this end we obtain analytical 
expressions for the lepton angular distribution, with 
a fixed value of the photon energy and general polar¬ 
ization. We then fold this expression with the photon 
luminosity function and the polarization profile for 
the ideal back-scattered laser spectrum 13 ; we obtain 
numerical results for the different mixed polarisation 
and charge asymmetries which we construct. Our choice 
of the ideal case for the back-scattered laser spectrum 
is for demonstration purposes. Further results using the 
recently available spectra, including the detector simu¬ 


lation for TESLA 


mectra, 

0 0 


will be presented elsewhere 


We then use the above-mentioned asymmetries 


P- y (s)-e IJjUa0 e t /e 2 k^kt/ 


( 2 ) 


where k\ and k 2 are the four-momenta of colliding 
photons and ei j2 are corresponding polarization vectors. 
We take S y , P 7 to be complex whereas S t ., Pt are taken to 
be real. This choice means that we assume only the CP 
mixing coming from the loop-induced effects in the Higgs 
potential. We allow these form factors to be slowly vary¬ 
ing functions of the 77 c.m. energy since in any model 
the loop-induced couplings will have such a dependence. 
Simultaneous non-zero values for P and S form factors 
signal CP violation. We will construct various asym¬ 
metries which can give information on these form factors. 

We allow the tbW vertex to be completely general 
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and write it as 

T * bW = ~7 5 Vtb [r UilPl + flRPR) 

- ^(Pi ~ Pb)u {f2 lPl + h rPr) , 
= -jf* b* (A lPl + IirPr) 

- - Pb)v (. ?2lPl + IirPr) , 


(3) 

(4) 


where Vtb is the CKM matrix element and g is the 
SU(2) coupling. We work in the approximation of 
vanishing b mass. Hence fm, fm, / 2 L and f- 2 R do not 
contribute. We choose SM values for f\L, and /il, 
viz., fiL = f\L = 1. The only non-standard part of 
the tbW vertex which gives non-zero contribution then 
corresponds to the terms with f 2 R and f 2 R. We will 
sometimes also use the notation f + and f~ for f 2 R 
and f 2 L respectively. One expects these unknown f’s 
to be small and we retain only linear terms in them 
while calculating the amplitudes. Below we give the 
helicity amplitudes for the production of tt followed by 
the decays of the t/t in terms of these general couplings. 


when both colliding photons have the same helicities, 
whereas the SM contribution is small for this combina¬ 
tion as we move away from the tt threshold. Thus a 
choice of equal helicities for both colliding photons can 
maximize polarization asymmetries for the produced tt 
pair, better reflecting the CP-violating nature of the 
s-channel contribution. It should be mentioned here 
that these statements are true only in the leading order 
(LO). Radiative corrections to 77 —> qq can be large 
0 0 That is also the reason we have restricted 
our analysis to asymmetries, which involve ratios. As 
a result the analysis is quite robust even if we use 
only the LO result for the SM contribution. Note also 
that the SM contribution for equal photon helicities is 
peaked in the forward and backward directions, whereas 
the scalar-exchange contribution is independent of the 
production angle 9 t . This also suggests that one can 
optimize the asymmetries by angular cuts to reduce 
the SM contributions to the integrated cross-section, 
of course taking care that the total event rate is not 
reduced too much. We will make use of this feature in 
our studies. 


A. Production Helicity Amplitude 


B. Decay Helicity Amplitudes 


The production process 77 —> tt receives the t/u chan¬ 
nel SM contribution from the first two diagrams of Fig. 1, 
which is CP-conserving whereas the s channel <f> exchange 
contribution may be potentially CP violating. The he¬ 
licity amplitudes for the s and the t/u channel diagrams 
are given by Eqs. (|[) and (||) respectively : 


We assume that the t quark decays only via the 
tbW vertex followed by the decay of W into lepton 
and corresponding neutrino. The helicity amplitudes 
Mr (At, A b, A 1 +, X„) and Mr(Af, A&, A;-, A p ), for the decay 
of t and t are given below: 


,, ,. , , , , —ieamt s 

M^Ai, A 2 ; A t , A t ) —-y—:- 

A-kMw s - + im^l $ 

[Sy{s) + *AiP 7 (s)] [A t(3S t - iP t ) ( 5 ) 

, . -iAnaQ 2 

Msm\Xi, X 2 ; X t , X t ) — - - -2 - 

1 — /P cos^ 9 t 

—7=- (Ai + At/?) 

. v s 

- -A A t /? sin 2 9 t Sx 1 ,-x 2 ^XtXt 

v s 

- 2/? (cos 0 t T AiAt) sin 9 t £x,,-A 2 £\ t ,-A r ] ■ (6) 


M r {+,-,+,-) = -i2V2g 2 y/m t plplp°+ A w (pw) 


, f2R.\ 9 l+ 

1 + 7?) ~ 


■ ®b i( h h 

cos — sin — e Vb — 


• 9 1 / 9 b j, 

sm — cos — e^ 1 ' 
2 2 

sin — sin -!A e i (? v -? > i+) + 


hR . Ob i 

-— sm —e 

^ 2 

0 V 9 t + 

cos — cos — 

2 2 


4 >b 


(7) 


Here /?, Q and 9t are velocity, electric charge and 
scattering angle of the t quark respectively; P^ and 
denote the total decay width and mass of the scalar 
</>; Ai j2 stand for helicities of two photons while the 
other A’s stand for helicities of particles indicated by the 
subscript. For photons, helicities are written in units of 
h while for spin- 1/2 fermions they are in units of H/2. 

From the expressions in Eqs. (|) and (||) it is 
clear that the </> exchange diagram contributes only 


-, +, -) = -x2\f2g 2 m t p 0 b plp° l+ A w (Pw) 


1 + 


/2 r\ ■ 0i+ _ 


sm-e 

2 




9 V . 9b 
cos — sm —e 
2 2 


i(j>b 


— sm — cos —e 
2 2 


i(/>n 


Hr Ob 

—=■ cos — 
2 


Ov . 0i+ i((h —d) ,) 0 U 9i+ 

sm — sm- e ’ + cos — cos- 

2 2 2 2 


( 8 ) 
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M r (+, +, +) = -i2\f2g 2 m t pffip°_ A w {p 2 w ) 


1 + 


f2L 


COS ■ 


up . v b __ i4> 
cos — sin —e Vb 
2 2 


— sin — cos — e 


$2l . 6b 


_■ sm — e z ^ b 
fr 2 


u v . u l~ -i(d>n—ch, ) I 6v 9i- 

— sin- p. vi-) \ ms — ms- 


sm — sm ——e 
2 2 


F ‘~ ’ + cos — cos ■ 

2 2 


, (9) 


M r {~, +, +) = -i2V2g 2 m t p^ b plp°_ A w {Pw) 


f2L\ . 6i 


sm — e l<t>l ~ 


9p . 9 b —i<ph 

cos — sm e Vb 


. 6* Ob -H,- 
— sm — cos — e 


/ 2 L 

—= COS — 

v'r 2 


sin -A sin -A-e l ^ v ^ l ~ ■* + cos -A cos 


0 t - 


,( 10 ) 


where, 


AwiPw) p 2^ _ Myy + iM w T w ’ r ml ' 

For simplicity, the above expressions for the decay ampli¬ 
tudes have been calculated in the rest frame of t (t) with 
the 2 -axis pointing in the direction of its momentum in 
the 77 c.m. frame. We treat the decay lepton l and the b 
quark as massless and list only the non-zero amplitudes. 


III. ANGULAR DISTRIBUTION OF LEPTONS 

Using the narrow-width approximation for the t quark 
and the W boson, the differential cross section for 77 —1 
ti —> l + bvit can be written in terms of the density matri¬ 
ces as 

da 3 e 4 g 4 (3Ei+ 

dcosdt d cos &i+ dEi+ d(j>i+ 64(47r) 4 sr t TO t r^M^ 

, + rr'(A,A')l 

P (AA) -^ 0 — ■ 

-- v -/ [ m tEf + J 

c.m. frame ^ v ^ 
rest frame 

(ii) 

Here E? + is the energy of l + in the rest frame of the 
t quark; the production and decay density matrices are 
given by 

P + (A, A') = e 4 p + (A, A') = ^ pi(Ai, A' 1 )p 2 (A 2 , A' 2 ) 
x M(A 1 ,A 2 ,A,A f )A4*(A' 1 ,A',A / ,A t -) 


E 

A A' 


T'(A,A') 

_ m t E° + 

rest frame 


and 

F(A, A') = g 4 \A(p 2 w )\ 2 T'(A,A') = E J da 

x ^ Mr(A, A b, A;+, A„) , A b , A;+, A„). 


Here a is the azimuthal angle of the b quark in the rest 
frame of the t quark with the 2 -axis pointing in the direc¬ 
tion of the momentum of the lepton. All repeated indices 
of matrix elements and density matrices are summed 
over; Pi( 2 ) are the photon density matrices; written, in 
terms of the Stokes parameters rg , : 


Pl(Al, A() : 
p 2 (A 2 , A 2 ) 


1 

2 


I + 72 

-T]3 - ii 71 


1 

2 


1 + 6 
-6 - *£ 1 


-m + ir) 1 
1 - 72 

-6 + *6 
1-6 


( 12 ) 

(13) 


Here, ry 2 is the degree of circular polarization while rft 
and 773 are degrees of linear polarizations in two trans¬ 
verse directions of one photon; 6 are similarly the degrees 
of polarization for the second photon. The explicit ex¬ 
pressions for the production density matrix p(\ tl \' t ) de¬ 
pend upon the polarization of initial photons. The decay 
density matrix elements are independent of any initial 
condition and in the rest frame of t (f) they are given by: 


r(±,±) 


r(±,T) 


r(±,±) 


r(±,T) 


g 4 m t E° l+ \A w {p 2 w )\ 2 (m 2 - 2p t .p t +) 

( 1 ± co; s 9 i+) (' 1 + * 1 A*)«L 
v y/r Pt-Pi+. 

g 4 m t Ei+\A w {pw)\ 2 {m 2 t - 2p t .p l+ ) 

sin d;+e ±i< ^ i+ (1 + m2R) M ™ ) , 
V v r Pt-Pi+J 


g 4 mtE | Aw (pw ) | 
(1 ± cos 61 -) ^1 + 
g 4 m. t E?-\A w (pw)\ 

sin Oi-e^i- (l + 


2 ( m 2 t - 2pi.pi-) 

m~2 l) 2 w A 

\fr Pi-Pi- J ’ 

2 (m 2 t - 2 pt.pi-) 

W2l) M 2 , A 
Vr Pi-Pi- J ' 


(14) 


(15) 


(16) 


(17) 


We have kept only the linear terms in the form factors 
f‘ 2 R and / 2 £, as we assume them to be small. These 
expressions written in terms of the lab variables can also 
be written in terms of the variables in the 77 c.m. frame. 
The relations between the angles in the rest frame and 
the 77 c.m. frame can be easily derived and are given by 


V 1 


sind;+e*^ i+ = — 

1 ... 

(sin Off 71 ' cos 0f' m ' cos rt+ 

+* sin Off™' sin< f>f +"'), 

sin 6 > ; -e^- = y/l ~^ 

1 - /3 cos 

(- sin Om cos 9f m - cos + cos 9 z c 

+* sin Ofd 71 ' smtfL™ ). 


j3 cos 9\™' 

ft? 1 - -cos 6 » 6 m ' sin Of 


Ofi m - s 



5 


Using the above relations and dropping the superscripts 
c.m. from the angles, we can rewrite Eq. © as 


da 


3a 4 /3 E l± 


d cos 9 t d cos 9i± dEi± d<j>i± 16i 2 frs T t T w M w 
1 4 Ei± 


8j t \ 1 - /? cos da frs{l~(3 2 ) 
3?(/±) 


x 1 + 




y/r E l± ^s(l - 13 cos 0 t i) 


x [^4 ± (1 - P cos 0 t i) ± B ± (cos0 t i - /?) 
tcVi -/? 2 sin &i± (cos 6 t cos </>;± — sin 9 t cot 9i± ) 


±D ± y/l^ 1 3 2 sin 0 ;± sin<fo±l , 


( 22 ) 


where = 1/-^/l — /3 2 , x w = sin 2 $vv, 9w being the 


Weinberg angle and 

A ± = p±(+, + )+p±(_,-), (23) 

B ± = p ± (+,+)-p ± (-,-), (24) 

C' ± = 2R[p *(+,-)], (25) 

D ± = -23[p' ± (+,-)]. (26) 


In Eq. (^) and what follows, the lepton variables are 
defined in the 77 c.m. frame. Further, 9ti stands for 
9 tl + for the l + distribution and hence the upper sign in 
the equation, whereas it stands for 9n~ for the lower 
sign and hence the l~ distribution. To get the angular 
distribution of leptons we still have to integrate Eq. ( |2^ ) 
over Ei, cos 9t and <j>i. The limits on Ei integration are, 


Myy 1 

yfs 1 - /3 cos 9 t i 


< Ei < -i 


1 


s 1 — /3 cos 9 t i 


After the Ei integration, we get 


and thus exactly cancels the one in Eq. (|2?]). Thus we 
see that the angular distribution of the decay lepton is 
unaltered, in the linear approximation of anomalous tbW 
couplings. In fact this is quite a general result, which is 
attained under certain assumptions and approximations 
we have made. We elaborate on this point below. 

An inspection of Eqs. (0M0) shows that the 
presence of any anomalous part in the tbW coupling 
changes the decay density matrix only by an overall 
energy-dependent factor independent of angle. The 
quantity Pt-Pi does have an apparent dependence on 
the angular variables of the lepton. However, in fact it 
depends only on the lepton energy. To see this clearly, 
let us go to the rest frame of the t quark. Now the 
three lepton variables are {E[ est , 0 [ est , (/>[ est } and the 
anomalous term depends only on £) rest . This means that 
the angular distribution of leptons in the rest frame of 
the t quark is unaltered by the presence of an anomalous 
term in tbW coupling, apart from an overall scaling. 
The angular distribution in any other frame can be 
obtained from that in the rest frame by a Lorentz boost. 
Thus the angular distribution of leptons in an arbitrary 
frame will be the same as that in the absence of the 
anomalous term, up to some overall factor that depends 
upon energy and the boost parameters and no angular 
variables. Of course, it is not very obvious by looking 
at Eq. (21) that this will indeed happen. But, with a 
change of variable, 

Ei -> E[ est = itEiil-P cos 0 tJ ), 
the additional overall factor becomes 

yfr m t Ej est ’ 


da 3ai 4 /3 1 

dcos9 t dcos9i± d<j>i± ldx^frs T t TwMw 
1 mf (1 + 2)—63?(/ ± ) v / r)(l — r ) 2 
87 1 6 s (1-/3 cos 9u) 3 


x [A ± (l - Pcos9 u ) ± B ± (cos9 t i - /3) 
+C ± V / 1©3 2 sin 9i± (cos 9 t cos </> ; ± — sin 9 t cot 9 [±) 


±D ± y/l- (3 2 sin 9[± sin(/);± 


(27) 


Here we have used the notation /+ and f~ for / 2 r and 
/ 2 L, respectively. From the above equation it is clear 
that the angular distribution of leptons after energy in¬ 
tegration is modified due to the anomalous tbW coupling 
only up to an overall factor 1 + 2r — 63?(/ ± ) v / r, which 
is independent of any kinematical variable. In fact, the 
same factor appears in the total width of the t quark 
calculated up to linear order in f ± : 


r *(t) = 


9 “3 

rag 

192x%, T\yMyy 


-r ) 2 


[l + 2r-65i(/ ± )V^ 
(28) 


which is clearly independent of angular variables. After 
integration over E[ est , in the limit M^ v /2m t < £) rest < 
TO t /2 we get back Eq. (|2?|). The important point is 
that in proving the result we did not make any reference 
to the production density matrix and hence the result 
is very general and applicable to any 2 —> 2 process for 
ti pair production provided the following conditions are 
fulfilled: 

• we use the narrow-width approximation for t and 
W, 

• b, l, v are taken to be massless, 

• the only decay mode of t is t —> bW —> blv, and 

• the anomalous tbW coupling is small enough that 
one can work to linear approximation in it. 

For the case of e + e _ —> ti followed by subsequent t/t 
decay, this was observed earlier Jl^, Q. It was proved 
recently by two groups independently; for a two-photon 
initial state by Ohkuma |2lJ| , for an arbitrary two-body 
initial state in p 2 | and further keeping to 5 non-zero 
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in |||]. These derivations use the method developed 
by Tsai and collaborators |24) for incorporating the 
production and decay of a massive spin-half particle. 
Our derivation makes use of helicity amplitudes and 
provides an independent verification of these results. 
The result is very crucial for our present work as we 
now have an observable where the only source of the 
CP-violating asymmetry will be the production process. 

Thus we can analyse the Higgs CP property easily, 
as long as the anomalous part of the tbW couplings, / ± 
is small and the quadratic term can be neglected. If / ± 
is not small then we have to keep the quadratic terms 
in Eqs. ©-®D and the decay density matrices to this 
order are then given by J 2 ?J : 


P(±, ±) = g 4 m t E?+ \/S.w{Vw )\ 2 i m t - 2 Pt-Pi+) 

m + ) M 2 , 


(1 ± cos 9 i+ ) 1 + 


Vr Pt-Pi+ 


\f' r (1 ~F cos 9 [+) I 1 


L w 


2 Pt-Pi+ 


+ 1 / 


+ | 2 _ 


cos e l+ 


2 r{p t .pi+) 2 

P(±,=F) = g 4 m t E°+\A w (pl)\ 2 (to? - 2 \p t -Pi+) 

, m + ) mi 


I / 


+ |2 


sin 6 »,+ 6 ^+ 

1 - 


m? + M'A 


w 


+ 


V~r Pt-Pi+ 

Ml 


2 Pt-Pi+ 2 r(p t .p l+ ) 2 J 


(29) 


(30) 


r(A, A') will be given by similar expressions. Thus if / ± 
are not small they can modify the angular dependence of 
the decay density matrix in the rest frame of the t quark 
and hence in any other frame. In that case it will not 
be trivial to use angular distributions to study the CP 
property of the production process. In this work we will 
assume to be small and will neglect the quadratic 
terms in Eqs. ( |29| ) and (fl0|). 

Making the above-mentioned four assumptions, which 
are very reasonable indeed, we now go on to calculate the 
final angular distribution by integrating Eq. (pTj) over 


cos 9 t and 4>i. We obtain for the angular distribution: 

da 37 ia 2 B . 

-T- a — = -2 [2^oo(^300 - flyhoi) 

d cos 0i± 8 S 7 ? L 

+2^01(^310 — Pyhn) + 2ALq 2 (/ 3 2o — Pyl 3 2i) 
+2^22(^322 — Pyl 3 23) + 2^43(1324 — Pyl 325) 

o± ^ 

±-|i{-2(/ 3 i 0 - Pyhu) + (1 - P 2 ) E x i W 


±^2{-2 (/ 320 - pyhn) + (1 - P 2 ) 


2=0 

2 


P 


±^{-2(7 32 i - Pyl 322 ) + (1 - P 2 ) Y ** W 


2=0 

3 


P 

B ± 


2=1 

5 


±—-p-{—2(7 323 - p y l 324 ) + (1 ^ P 2 )J2 X * « 


P 

6 

±Co Y Y i J 52i 
7—0 


i- 3 


(31) 


where 


7 ijk — 


d cos 6 t cos k d t 


-1 


[a + (y — P cos 0 t) 2 ]®/ 2 (l — p 2 cos 2 0 t ) J 


a = (1 — P 2 ) sin 2 9i ±, 

V = cosd[±. 

We have obtained explicit analytical expressions 
for all Iijk appearing in Eq. ([fl]), which are not listed 
here. Aij and Bij are coefficients in expansions of the 
following type: 

„± v- A tj c os l 0t 
A =E(T77w7)J ^ 

*>J 

Expressions for d^’s and Bj 3 ’s for circular polarization 
of photons and expressions for XP s, Y)’s and Co are given 
in the appendix. Equation © is the angular distribu¬ 
tion of leptons for a given 77 centre-of-mass energy. In 
a 77 -collider constructed using the back-scattered laser 
beam one will not have monoenergetic photons in the 
inital state; further, the degree of circular polarization 
of the photons will depend on its energy. Thus the fi¬ 
nal observable cross section is to be obtained by folding 
Eq. (|3^) with the luminosity function after accounting 
for the energy dependence of the circular polarization of 
the photons. 


IV. PHOTON COLLIDER 

In a 77 collider, high energy photons are produced by 
Compton back-scattering of a laser from high energy e~ 
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or e + beam via 


e (A e -) 7 (A ; J -> e 7 (Ai) 
e + (A e +) 7(A; 2 ) -> e + 7(A 2 ). 


In this paper we will be using the ideal photon spectrum 
due to Ginzburg et al. for x < 4.8. The ideal luminosity 
(for zero conversion distance) is given by 


1 (Hj ^ ry 

L e - e + dyidy 2 


where 


f(y) = 


2-Kor 


1 


a c xm1 |_ 1 — j/ 

- \ e \irx(2r — 1)(2 — y)], 


= /(j/i)/( 2 / 2 ), 


+ 1 — y — 4r(l — r) 


AEb “° - 15.3 I 

( E b \ 

(-) 

m e 

\Ee\) 

\eVJ 


x(l - y) 


< 1 , 


V = W 

<J C = cr” P + A e A/CTi, 


U) < 


xE b 
1 + x' 


2ira 2 


xmi 


x x- 

1 


1 -o ) 1 °S ( 1 + X ) 


+ 2 + x 2(1 +cc) 2 J ’ 


27TOT 


O' 1 = 


1 + - log(l + x) 

X 


xm% 

5 1 1 

2 + 1 + x ~ 2(1 +a;) 2 


(32) 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


A e and A; are the initial electron (positron) and laser 
helicities respectively, and uiq is the energy of the laser. 
In Eq. (f32|), if we change variables from y\ and y 2 to 
z = ^Jy\yn and y 2 and integrate over y 2 , we will get an 
expression for the photon spectrum as a function of the 
77 invariant mass W = 2zE b , where E b is the energy 
of the e ± beam and is plotted in Fig. [ 2 ] for x = 4.8. 
The spectrum is peaked in the hard photon region for 
A e A; = —1. 


The mean helicity of high energy photons depends 
upon their energy in the lab frame. In an ideal 77 
collider the energy dependence of the mean helicity is 
given by 


27TO' 2 

Viy) = ^mg/(„) 1A ‘ (1 ~ 2 ’' )(1 ~ i/ + 1/(1 ~ i,)) 

+ A e ra[l + (1 - y)(l - 2r) 2 ]} (40) 

and is plotted in Fig. || for x = 4.8. For A e A; = — 1 
the back-scattered photon has the same helicity as the 
electron (positron). Also, the spectrum is peaked at high 
energy, and yields a high degree of polarisation of the 



FIG. 2: Luminosity distribution plotted against a (which is re¬ 
lated to the 77 invariant mass W = 2 - v /wiu >2 by z = W/(2E b )) 
for x = 4.8. Solid line corresponds to A e A; = —1, small dashed 
line is for A e Ai = 1 and large dashed line is for A e A; = 0. Con¬ 
version distance is taken to be zero. 



0.2 0.4 0.6 0.8 1 

aj/E b 


FIG. 3: Mean helicity of scattered high energy photon plotted 
against reduced energy of photon y = a j/E b . Solid lines are 
for A e Ai = —1 and the dashed lines corresponds to A e Ai = 1. 
The lines are marked with the values of (A e , Ai). 


photon beam. Hence, the dominant photon polarization 
in this case is decided by the electron (positron) helicity. 
Now, as suggested by Eq. (||), the helicities of two 
colliding photons should be equal in order to have 
Higgs contribution. Thus we choose A e A/ = —1 to get a 
hard photon spectrum, and set A e - = A e + to maximize 
the Higgs contribution, and hence the sensitivity to 




























possible CP-violating interactions. For our numerical 
analysis, we have chosen A e A; = —1; the initial state 
can thus be completely described by the helicities of the 
initial electron and positron. We denote the total cross 
section in the lab frame by <r(A e -,±), where the sec¬ 
ond argument denotes the charge of the final state lepton. 

The total cross-section with an angular cut in the 
lab frame can be obtained by folding Ecp © with the 
photon spectrum: 

g(6o,yi,V2) 

c'(A e -,±)= l f L y dyidm f d cos 0\± ^ , 

J dyidy 2 J dcosO[± 

f(6o,yi,y2) 

(41) 

where f(0o, yi,y 2 ) and g(Oo, 2 /i) 2 / 2 ) are the (boosted) lim¬ 
its on integration in the 77 c.m. frame. We end this 
section with a few remarks. We have presented the spe¬ 
cific case where the 77 collider is based on a parent e + e - 
collider and we also assume 100 % polarization for the 
e _ /e + . The analysis is completely valid for the case of a 
parent e _ e _ collider, for which achieving a high degree 
of polarization for initial leptons might be technologically 
simpler. 

V. NUMERICAL RESULTS 

To determine the CP properties of the Higgs, we need 
to know all the four form factors appearing in Eqs. (l|) 
and (|[). Assuming the mass and decay width of Higgs 
to be known, we then have the following six unknowns 

S u Pt, »(S 7 ), 9f(S 7 ), 3?(P 7 ), 3(P 7 ). 

They appear in eight combinations in the expression for 
the production density matrix, which we denote by Xi 
and y-i (i = 1,4), and are listed below, together with 
their CP properties. 


Combinations Aliases CP property 


St^(S 7 ) 

Xi 

even 

St$*(S 7 ) 

X 2 

even 

St*t{P-y) 

2/i 

odd 

smpj 

2/2 

odd 

PMS-y) 

2/3 

odd 

PMS-y) 

2/4 

odd 

P t 3?(P 7 ) 

x 3 

even 

PMP-r) 

X\ 

even 


Only five of these combinations are independent because 
they satisfy the following relations 

2 /i ■ Vs = Xi ■ x 3 , 2/2 ■ 2/4 = x 2 . £4, 

yi . x 4 = 2/2 • x 3 , 2/4 • xi = 2/3 ■ x 2 . 

Any three relations of the four listed above are in¬ 
dependent relations while the fourth one is derived. 


Expressions for asymmetries can by written in terms of 
x’s and y’s and can be used to put limits on sizes on 
these combinations. 

In what follows we will define various asymmetries 
involving the polarization of the initial e and charge of 
the final decay lepton, some of which are CP-violating 
and use them to put limits on the size of various 
combinations of the form factors. There is no forward- 
backward asymmetry because two photons with the 
same helicities are indistinguishable in their c.m. frame. 
That is, no favoured direction exists and the forward 
direction is indistinguishable from the backward. This 
is to be contrasted with the situation studied in [p5| , 
where forward-backward asymmetry could be used to 
put limits on CP violation arising from the top electric 
dipole moment or a CP-odd 77 Z coupling. The effects 
of the s-channel Higgs-exchange diagram appear only in 
charge and polarization asymmetries along with purely 
CP-violating asymmetries. For our numerical studies 
we take the values of the form factors calculated in the 
MSSM for certain values of its parameters. The specific 
values which we use for demonstration purposes are 
taken from Ref. 0. These are for tan/3 = 3, with all 
sparticles heavy and maximal phase: 

= 500 GeV , = 1.9 GeV, 

S t = 0.33 , P t = 0.15, (42) 

S 1 = -1.3-1.2 i , P 7 = -0.51 + l.U. 

For the SM, 5”s and P’s are identically zero. By SM 

we mean contribution only from t and u channels. The 
light CP-even Higgs contribution at the tt threshold and 
beyond is small and hence is neglected. 

A. Polarized Cross Sections and Asymmetries 

There are two possibilities for the initial state polar¬ 
ization, A e - = A e + = +1 and —1. In the final state we 
can look for either l + or l~. This makes four possible 
polarized cross sections listed as 

n(+,+), cr(+, — ), er(—,+), cr(—, —). 

These are plotted in Fig. |as a function of the electron 
beam energy Eb for the angular cut of 60° in the lab 
frame. For the SM, cr(+,+) and er(—,—) are exactly 
equal as they are CP conjugates of each other. In 
the MSSM, because of CP violation, they can differ. 
A similar statement can be made about the pair 
(t(—, +), er(+, —). The flat behaviour with energy of 
curves 3 and 4 is due to the destructive interference 
of the Higgs-mediated amplitude with the continuum. 
Recall here again that second index in the expressions of 
the cross sections is the sign of the charge of the lepton. 
A comparison of curves 1,3 and 5, then shows clearly 
the change in the sign of interference effects as the sign 
of polarizations of the two photons is changed from ++ 
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FIG. 4: All four integrated cross sections are plotted against 
the beam electron energy Eb , for the SM as well as for MSSM 
with our choice of parameters. The angular cut used in this 
figure is 60° in the lab frame. Line 1 is for <r(+,+), line 
2 for er(+, —), line 3 for <r(—, —) and line 4 for er(—,+) in 
the MSSM. Line 5 is for <r(+, +) and <r(—,—) and line 6 for 
ct(+, —) and er(—, +) in the SM. 


to-. The jump in ct(+, +) and ct(+, —) at around 

310 GeV corresponds to matching of Higgs resonance 
peak with the peak of the hard photon spectrum. This 
suggested to us the choice of E b = 310 GeV for the 
analysis, as the deviation from the SM is then very large 
for the chosen value of parameters. 


We choose two polarized cross sections at a time, 
out of the four available, and define six asymmetries as 


a _ a (+> +) - cr(—, ~) 
ct(+,+) + ct(-,-)’ 

, _ Q'C+i ~) ~ &(—> +) 

2 ct ( + ,-)+ ct (-,+)’ 

A _ g ( + i +) ~ g (~1 +) 

cr(+, +) + o(—, +) ’ 
a = tj(+,~) - cr(-,-) 
ct(+,-) + ct(-,-)’ 

A _ g ( + ) +) ~ g ( + i ~) 

5 ct (+,+) + ct (+,-)’ 

A = g(~,+) - 0-(—,-) 

6 a(-,+) + a(-,-)- 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


of the lepton, and As and Aq are charge asymmetries for 
a given polarization. All these asymmetries are plotted 
against the beam energy E b for SM and MSSM in Fig. [|. 
From these plots it is clear that E b = 310 GeV is a good 
choice for putting limits on the size of the form factor, 
for our choice of the mass of the scalar . 


B. Sensitivity and Limits 


After choosing a suitable beam energy for the analysis, 
the next thing to look for is a suitable angular cut in the 
lab frame, which will maximise the sensitivity of the mea¬ 
surement. For asymmetries to be observable, the number 
of events corresponding to the asymmetry must be larger 
than the statistical fluctuation in the measurement of the 
total number of events. If N is the total number of events 
then the number of events corresponding to the asymme¬ 
try must be at least f\/N, where / = 1.96 for 95% C.L. 
The number of events N = o £, where C is the luminosity. 
Asymmetries are defined as 


A- a 1 ~ 

CTl + 0 2 


Act 

ct 


Thus the number of events corresponding to the asym¬ 
metry is CAo. For the asymmetry to be measurable at 
all we must have at least CAo > fVCo , with / de¬ 
noting the degree of significance with which we could 
assert the existence of an asymmetry. Thus the ratio 
(CAo)/{fV~Co) = {VC/f) x (A o/y/o) will be a measure 
of the sensitivity. One can be more precise in defining 
this by noting that the fluctuation in the asymmetry is 
given by 

for A <C 1. The larger the asymmetry with respect to the 
fluctuations, the larger will be the sensitivity with which 
it can be measured. We define sensitivity as, 

c A Act 

5 = 77 * 

5A \fo 


Ao/a/o, which is proportional to the sensitivity, is 
plotted for all asymmetries in Fig. |(i| against the angular 
cut in the lab frame, 9q. Since A\ and A 2 are purely 
CP-violating, they have no contribution from SM for 
any angular cut. Hence for them, the sensitivity is 
large when the angular cut is small, because of better 
statistics. For the other four there is an SM contribution 
that varies with the cut. Though the exact position of 
the peak in Fig. depends upon the relative sizes and 
signs of the form factors, 9q = 60° seems to be a good 
choice for the angular cut to maximise the sensitivity of 
four of the asymmetries. 


Of the above six, A\ and A 2 are purely CP-violating, A 3 The process under consideration violates CP in general, 
and A 4 are polarization asymmetries for a given charge But when the cut 9g is 0, the partial cross sections for 
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FIG. 5: All six asymmetries are plotted as a function of beam energy Eb for the SM (dashed line) and the MSSM (solid line) 
at an angular cut of 0o = 60° in the lab frame. At Eb = 310 GeV, owing to resonance in the s-channel, the MSSM values of 
asymmetries are maximally different from that of the SM. 
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FIG. 6: Aa/y 7 ^, which is proportional to the sensitivity S, is plotted against the angular cut do for all six asymmetries. Solid 
line is for the MSSM and dashed line for the SM. For Mi and M 2 , the smaller angular cut is favourable while for others do = 60° 
is a good choice. 


/ + and l~ production become the corresponding total 
cross sections, and are therefore equal, because of charge 
conservation. Hence for 8 q —> 0, A 5 and Me approach 
zero. In that limit, the polarization asymmetries M 3 
and A 4 are purely CP-violating. Thus for M 3 and M 4 , 


apart from the choice 9q = 60°, where the sensitivity 
peaks, 0 o = 0 ° would also be a good choice for isolating 
CP-violating parameters. But, to be away from the 
beam pipe, we choose the lowest cut to be 20 ° in the 
lab frame. We have used all six asymmetries for angular 
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TABLE I: List of 95% C.L. limits on all the combinations at 
500 fb" 1 and 1000 fb" 1 at Eb = 310 GeV. These limits are 
obtained from data plotted in Fig. [t| 



min 

max 

min 

max 

MSSM 


(500 fb -1 ) 

(500 fb -1 ) 

(1000 fb -1 ) 

(1000 fb" 1 ) 

value 

Xi 

-3.775 

3.594 

-2.990 

2.869 

-0.429 

X2 

-3.413 

3.896 

-2.748 

3.111 

-0.396 

X3 

-2.386 

2.873 

-1.842 

2.386 

-0.077 

#4 

-2.837 

2.465 

-2.375 

1.930 

+0.165 

2/i 

-2.786 

2.786 

-2.148 

2.148 

-0.168 

2/2 

-3.095 

3.095 

-2.433 

2.433 

+0.363 

2/3 

-2.155 

2.155 

-1.687 

1.687 

-0.195 

2/4 

-2.346 

2.346 

-1.867 

1.867 

-0.180 


cuts of 20 ° and 60° to put limits on the combinations 
Xi and yi. If for certain values of the form factors the 
asymmetries lie within the fluctuation from their SM 
values, then that particular point in the parameter space 
cannot be distinguished from SM at a given luminosity. 
That point will be said to fall in the blind region of the 
parameter space. Thus the set of parameters {xi,yi} 
will be inside the blind region at a given luminosity if 

\-A({xi, //»})—AIsm| < £4sm = ^J sm £ V^ + ^m- ( 49 ) 

For simplicity we have taken only two of the eight com¬ 
binations to be non-zero at a time and have constrained 
them in 16 different planes, shown in Fig. |t], satisfying 
their inter-relations. The limits obtained on each of the 
combinations by taking a union of the blind regions in 
the 16 plots are listed in Table ||. Also shown in the 
last column of the table are the values of Xj'y 3 for the 
MSSM point we have chosen for Figs. || and |(]. 

Next we do a small exercise to see whether these 
asymmetries have the potential to distinguish between 
the SM and MSSM. It is clear that we can repeat the 
analysis of finding blind regions in the ( x 3 ,yj ) planes 
around a particular point predicted by the MSSM. 
The values of Xi,yj, corresponding to our choice of the 
MSSM point given by Eq. (^) are listed in last column 
of Table jjj. The blind regions around these values will be 
defined by an equation similar to Eq. (^), where AIsm 
will be replaced by AImssm- In A(xk,yi) all independent 
Xk, yi other than the pair a+, yj being considered, are set 
to their MSSM value, and the pair Xi,y 3 is then varied. 
We show in Fig. || the results of such an exercise for the 
parameter pair (£ 3 , 2 / 3 ) along with the corresponding 
one for the SM. This shows that these studies can be 
sensitive to the CP mixing produced by loop effects. Of 
course one needs to study this over the supersymmetric 
parameter space. But the example shown here clearly 
shows the promise of the method. 


VI. DISCUSSION 

The four cross sections depending on the polarization 
of the initial lepton and the charge of the final state lep¬ 
ton that we use to construct asymmetries, can in general 
be written as 

cr{X e ,Qi) = ctoo + aoiQi + cr 10 A e + anQiX e . (50) 

This says that we have four independent aij, which 
constitute four polarized cross-sections. Out of these 
four, fj 00 is the largest and others are of the order of a 
few per cent of croo- Thus we can safely approximate 
denominators of all A, 's to a = 2cr 0 o • This makes Aids 
proportional to their numerators, which consists of only 
three of aij. Thus out of six asymmetries constructed 
in Section [V] only three are independent and we cannot 
determine all six form factors simultaneously using these 
asymmetries. This is a reflection of the fact that there 
are only three CP-violating asymmetries |2(| at the 
production level of the tt pair; one is for the unpolarized 
case, and the other two are polarization asymmetries. 
The Ads defined here are combinations of these three. 

In Fig. [i] we took only two combinations as non¬ 
zero and varied them to find the blind region in that 
plane. We found strong limits on y’s and almost no 
limits on ads in each of the planes. When we allowed 
three of the combinations to vary simultaneously there 
were almost no limits on any of the combinations. This 
can be understood by looking at Eq. (|5(]). The charge 
asymmetries are very small and approach zero as we 
reduce the angular cut, which implies that 001 and an 
are very small and tend to zero as $o —> 0- Thus two 
of the four independent components of the polarized 
cross-section are very small (typically by a factor of 
100 to 500); neglecting them, we are left with only two 
independent components, implying that only one of the 
six asymmetries is independent. Thus, though we have 
four independent components at hand, two being small 
we are effectively left with only two almost identical 
strong constraints, and thus essentially only one. These 
asymmetries thus constrain only y’s and leave ads mostly 
unconstrained. The fact that ads are constrained at all is 
because the equation of the boundary of a blind region 
arising from any one asymmetry, for two variables, 
is an equation of a pair of conic sections. The blind 
regions shown in Fig. |?] are intersections of blind regions 
obtained from all six asymmetries with two different 
angular cuts. 


A. The Strategy 

All the cross sections and asymmetries are expressible 
in terms of ads and y’s alone and, of these, only five are 
independent. Thus any number of asymmetries for any 
general polarization can never determine all six form 
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FIG. 7: The boundaries of blind regions in the parameter space are plotted for various pairs of parameters, for luminosities 
500 and 1000 fb -1 at beam energy Eb = 310 GeV. Both angular cuts, 9q = 20° and 60°, are used to put limits at C.L. of 95%. 
The larger region corresponds to 500 fb 1 , while the smaller corresponds to 1000 fb -1 . 



FIG. 8: The boundaries of blind regions in the parameter 
space are plotted in xz~yz plane, for a luminosity of 1000 fb -1 
at beam energy Eb = 310 GeV. Both angular cuts, do = 20° 
and 60°, are used to put limits at C.L. of 95%. For this MSSM 
point, xz = —0.077, yz = —0.195. 


factors as only five independent combinations appear in 
the expressions. For St and Pt we have to rely on partial 
decay width measurements of the scalar <f> to tt pair. 
Thus, if we have a few more independent and strong 
constraints, we will be able to put simultaneous limits 
on all six form factors. But with circular polarization 
we have only the four observables used here already. 
One possibility would be to use the dependence of the 
angular distribution of the decay leptons on the initial 
state photon polarization. But to do that one would 
need a large statistics, which will not be available even 
with an integrated luminosity of 10 3 fb -1 . The other 
option is to look for linearly polarized initial photons. 
Here, by choosing different angles between the planes 
of polarizations one can alter the relative contribution 
from CP-even and CP-odd Higgs. This, along with the 
asymmetries considered and the partial decay width of 
cA, can then be used to put limits on all six form factors 
simultaneously or alternatively to determine them. Some 
discussions of these for the tt production exist already j?j] 
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FIG. 9: Blind regions along the line joining two model points, 
SM and MSSM, in the five-dimensional parameter space at C 
= 500, 750 and 1000 fb” 1 . All six asymmetries with both 
angular cuts, do = 20° and 60°, are used. 


In view of above analysis we can propose a strat¬ 
egy for characterizing the heavy scalar <f>. The first step 
would be to determine its mass m^, its total decay width 
and the partial decay width to a tt pair. The last will 
tell us about <Sf + P t 2 . Then the second step will be to 
look for asymmetries Ai and A 2 to see if there is any CP 
violation. Step 3 depends upon the outcome of step 2. 
In case of non-observation of CP violation, one will have 
to look for linearly polarized asymmetries to see whether 
the Higgs is CP-even or CP-odd. If CP violation is 
observed, then all the asymmetries, for circular and 
linear polarizations, can be used to determine the form 
factors. 


B. Discriminating Models 

As we have seen, it is not possible to determine all 
the combinations Xi , yj using the asymmetries we have 
constructed. However, as discussed below, we can surely 
use the model predictions for these to discriminate 
against a particular model when data are available or 
test the possibilities of being able to distinguish between 
different models at a given luminosity. 

The blind region around any model point in the 
five-dimensional parameter space is a non-convex struc¬ 
ture and extends far out from the model point in some 
of the directions. Thus projection on any plane may 
result in a large blind region, which can be misleading. 
Thus it is not possible to restrict to less than the full set 
of 5 parameters for testing models. Below we develop a 
method for distinguishing between models and checking 
whether they are ruled out by experiment. 

The simplest way to compare two models is to ask 
if the first model point lies within the blind region 
of the second and vice versa. If not, we say that the 
model predictions are distinguishable from each other 
at the luminosity and confidence level considered. As 
an example, we chose two models - SM and MSSM. The 
MSSM is same as given by the last column of Table 
|jj The model points in the five-dimensional parameter 


TABLE II: The probability V of confusing SM with MSSM 
at 95% C.L. for different asymmetries and luminosities. 


Asymmetries V at 500 fb 1 P at 750 fb 1 V at 1000 fb 1 


Ai{ 9 q 

= 60 u ) 

5.5 x 10 ” 4 

4.8 x 10 “° 

6 x 10” 8 

^2(^0 

= 60 °) 

7.8 x 1 CT 4 

8.0 x 10” 6 

7 x 10 ” 8 

A3 (6 0 

= 60 °) 

3.4 x 1 CT 4 

2.0 x 10” 6 

1 x 10” 8 

^4(^0 

= 60 °) 

1.4 x 1 CT 3 

2.1 x 10” 6 

2.8 x 10” 7 

Ai(6*o 

= 20°) 

1.6 x 1 CT 7 

< 10” 8 

< 10” 8 

^2(^0 

= 20°) 

2.3 x 10 ” 7 

< 10” 8 

< 10” 8 

A3 (6 0 

= 20°) 

2.3 x 10 ” 7 

< 10” 8 

< 10” 8 

A^Oq 

= 20°) 

2.2 x 10” 7 

< 10” 8 

< 10” 8 


space are connected by a line parameterized by t with 
t = 0 corresponding to one model and t = 1 to the 
other. We have calculated the blind region around each 
of the models along the connecting line. These are 
shown in Fig. ^ and it can be seen that each model 
sits well outside the blind region of the other at an 
integrated luminosity of 500 fb -1 . Furthermore, their 
blind regions do not overlap along these lines. Thus 
we can say that our method can distinguish candidate 
models at a certain luminosity (500 fb -1 in this case). 
A more accurate way will be to search the whole of the 
five-dimensional parameter space for the overlap of the 
blind regions corresponding to two candidate models 
and not just along the line joining them. If no such 
overlap is found then we can say for sure that the models 
can be distinguished. This search could be quite com¬ 
plicated. Alternatively we can use the numerical values 
of individual asymmetries and fluctuations directly as 
discussed below. 

It is clear that we can determine blind region around a 
given model prediction in any parameter space given the 
numerical value of the model predictions for asymmetries 
and the statistical fluctuations expected in it at a given 
luminosity. Any change in these numerical values will 
yield a different blind region in the five-dimensional 
parameter space. One will then test asymmetry pre¬ 
dictions for a particular model against an experimental 
measurement or compare the predictions of two models 
against each other to draw conclusions about their 
distinguishability at a given luminosity and confidence 
level. If the values of asymmetries expected at the 
particular level of confidence, corresponding to (say) 
two models, have no overlap, then the two models are 
distinguishable at that confidence level. There is still a 
non-zero probability that the models can be confused 
with each other in an experiment. To determine the 
probability of such a confusion we take any one asym¬ 
metry at a time and calculate the limits upto which the 
predicted asymmetry values can fluctuate at a certain 
level of confidence in the models under consideration. 
Then we generate normally distributed random numbers 
centered at the asymmetry corresponding to the first 
model, say SM, with standard deviation same as the 1 a 
















fluctuation of the SM asymmetry. We count the number 
of points for which the asymmetry value lies within the 
95% confidence interval of the other model, say MSSM. 
The number of such points divided by the total number 
of points taken is the probability V of confusing SM with 
MSSM at 95% confidence level. As Table [id| indicates, 

V is of the order of 10 -7 for a luminosity of 500 fb” 1 , 
and we can safely say that SM is distinguishable from 
MSSM at 500 fb -1 . In a similar way we can replace 
SM by the experimental asymmetries and MSSM by a 
candidate model. Now even if for one of the asymmetries 

V is very small (0(1CU 3 )) we can simply reject the 
model as in words it translates to: the probability of the 
experimental results being statistical fluctuation of the 
candidate model at 95% C.L. is very small. In fact, the 
method described above is nothing but the Step 2 of our 
strategy discussed in previous sections, when one talks 
about asymmetries Ai and ,4 2 . 
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VII. CONCLUSION 

We have demonstrated how the Higgs-mediated CP 
violation in the process 77 —> tt can be studied by look¬ 
ing at the integrated cross section of l + /l~ coming from 
the decay of t/t. We demonstrated that the decay lepton 
angular distribution is insensitive to any anomalous part 
of the tbW coupling, f ± , to first order. We constructed 
combined asymmetries involving the initial lepton (and 
hence the photon) polarization and the decay lepton 
charge. We showed that using only circularly polarized 
photons will be inadequate to determine or constrain 
the sizes of all form factors simultaneously, but can put 
strong limits on CP-violating combinations, y’s, when 
only two combinations are varied at a time. We show, 
by taking an example of a particular choice of MSSM 
parameters, that the analysis is sensitive to the CP 
mixing at a level that is generated by loop effects. We 
also further sketch a possible strategy to characterize 
the scalar <j> using linear polarization. 
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APPENDIX A: EXPRESSIONS FOR A ZJ , B , :j , ETC. 


For circularly polarized photons the form factors Aij 
and Bij are given below; 92 and £ 2 are the degrees of 
circular polarization of two colliding photons: 

A± = 2|A 0 | 2 (/3 2 5 t 2 + P t 2 )[(|5 7 | 2 + |P 7 | 2 ) 
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X 0 = 2 + /3 2 sin 2 61 
X\ = —4/3 cos 61 
X 2 = /3 2 (3cos 2 9i — 1) 

To = — cos 9i (2 + /3 2 sin 2 61 ) 

Y\ = /3(3 + cos 2 9i) 

Y 2 = 2 cos 9 1 (2 — (3 2 cos 2 9i) 

Y 3 = -2/3(3 + cos 2 9{) 

1 4 = cos 9i(—2 + 3/3 2 + /3 2 cos 2 9i) 
Y 5 = /3(3 + cos 2 9i) 

Yq = — 2/3 2 cos 9 1 


C 0 = 
D ± = 
A-c — 

Afj, = 


























15 


[1] J. A. Aguilar-Saavedra et al. fE CFA/DESY LC 
Physics Working Group Collaboration], |icp-pli/0106315| ; 
B. Badelek et al. [ECF A/DESY Photon Collider Working 
Group Collaboration], hep-ex/0108012 . 

[2] S. Weinberg, Phys. Rev. Lett. 37 (1976) 657. 

[3] T. D. Lee, Phys. Rep. C 9 (1974) 143. 

[4] A. Pilaft sis and C. E. M . Wagner, Nucl. Phys. B553 
(1999) 3 ( hep-ph/9902371 ); S. Y. Choi, M. Drees and J. 
S. Lee, Phys. Lett. B481 (20001 57; M. Carena, J. El¬ 


[15] V. I. Telnov, Nucl. Instrum. Meth. A 355 (1995) 3; A 
code for the simulation of luminosities and QED back¬ 
grounds at photon colliders, talk presented at Second 
Workshop of ECFA-DESY study, Saint-Malo, France, 
April 2002. 

[16] R.M. Godbole, S.D. Rindani and R.K. Singh, In prepa¬ 
ration. 

[17] G. Jikia and A. Tkabladze, Phys. Rev. D 54 (1996) 2030 
( i hep-ph/960 1384T); Phys. Rev. D 63 (2001) 074502 (ihep^ 


lis, A. Pilaftsis and C. E. M Wagner, Nucl. Phys. B586 
(20001 92 (Ihep-ph /0003 1800. S .Y. Choi and J .S. Lee. 
Phys. Rev. D62 (2000) 036005 ( |hcp-ph/991233c| ). 

[5] E. Asakawa, J.-i. Kamoshita, A. Suga moto and I. Wata n- 
abe, Eur. Phys. J. C 14 (20 00) 335 (|hcp-ph/9912373| ). 

[6] S. Bae, B. Chung and P. Ko, hep-ph/0205212; S.Y. Choi, 
B. Chung, P. Ko and J.S. Lee, [hcp-ph/0206025 . 

[7] E. Asakawa, S. Y. Choi, K . Hagiwara and J . S. Lee, Phys. 
Rev. D62 (2000) 115005 ( |hep-ph/0005313| ). 

[8] B. Grzadkowski and J. F. Gunion, Phys. Lett. B 294 

(19921 361-368. _ 


ph/0004068)7 


[18] M. Melles, W. J. S tirling and V. A. Khoze, Phys. Rev. D 
61 (2000) 054015 ( |hep-ph/990723^ ). 

[19] S.D. Rindani, Pramana 54 (2000) 791 ( hep-ph/0002006|). 

[20] B. Grzadkowski and Z. Hioki, Phys. Lett. B 476 (2000) 


87 (hep-ph/9911505); Z. Hioki, hep-ph/0104105. 

[21] K. Ohkuma. Nu cl.Phvs.Proc.Suvvl. Ill (2002) 285, 
(|hep-ph/0202126|). 


[9] G. Jikia and S. Soldner-Rembol d, in Ref. 1; Nucl. In- 
strum. Meth. A 472 (2001) 133 Qhep-ex/010105q ). 

[10] P. Niezurawski, A.F. Zarnecki and M. Krawczyk, hep- 
|ph/0207294|. 

[11] S. Y Choi, D. J. Mille r, M. M. Muhlleitner and P. M. Zer- 
was, bep-ph/0210077. 

1. AL 


[22] B. Grzadkowski and Z. Hioki, Phys. Lett. B 529 (2002) 
82, ( |hcp-ph/011236l| ). 

[231 B. Grzadk owski and Z. Hioki, F T-19-02. ( frrep- 


ph/0208079); Z. Hioki, hep-ph/0210221. 


[12] M. M. Muhlleitner, M. Kramer, M. Spira and P. M. Z er- 


was, Phys. Lett. B 508 (2001) 311 (Jiep-ph/0101083). 

[13] I. F. Ginzburg, G. L. Kotkin, S. L. Panfii, V. G. Serbo 
and V. I. Telnov, Nucl. Instrum. Meth. 294 (1984) 5. 

[14] A.F. Zarnecki, CompAZ: parametrization of the 


ICHEP’ 2002. abstract #156: 

hep-ex /0207021I 

http://info.fuw.edu.pl/~zarnecki/comp8 

iz/compaz.html[ 


[24] S.Y. Tsai, Phys. Rev. D 4 (1971) 2821; S. Kawasaki, 
T. Shirafuji and S.Y. Tsai, Prog. Theo. Phys. 49 (1973) 
1656. 

[25] P. Poulose and S.D. Rindani, Phys. Rev. D 57 (1998) 
5444, D 61 (2000) 119902 (E); Phys. Lett. B 452 (1999) 
347; P. Poulose, Nucl. Instrum. Meth. A 472 (2001) 195. 

[26] W.-G. Ma, C.-H Chang, X.-Q. Li, Z.-H. Yu and L. Han, 
Commun. Theor. Phys. 26 (1996) 455; Commun. Theor. 
Phys. 27 (1997) 101. 

[27] Expressions for T are written in the rest frame of the t 
quark. 



























































